Bulk Termination of the Quasicrystalline Five-Fold Surface of Al7oPd2iMn9 
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The structure of the Al7oPd2iMn9 surface has been investigated using high resolution scanning 
tunnelling microscopy (STM). From two large five- fold terraces on the surface in a short decorated 
Fibonacci sequence, atomically resolved surface images have been obtained. One of these terraces 
carries a rare local configuration in a form of a ring. The location of the corresponding sequence of 
terminations in the bulk model M of icosahedral i-AlPdMn based on the three-dimensional tiling 
q-*(2F) F- phase has been estimated using this ring configuration and the requirement from the 

LEED work of Gierer et al. that the average atomic density of the terminations is 0.136 atoms per 
A^. A termination contains two atomic plane layers separated by a vertical distance of 0.48A. The 
position of the bulk terminations is fixed within the layers of Bergman polytopes in the model M: 
they are 4.08 A in the direction of the bulk from a surface of the most dense Bergman layers. From 
the coding windows of the top planes in terminations in M we conclude that a Penrose (PI) tiling 
is possible on almost all five-fold terraces. The shortest edge of the tiling PI, is either 4.8A or 7.8A. 
The experimentally derived tiling of the surface with the ring configuration has an edge-length of 
8.0 ± 0.3 A and hence matches the minimal edge-length expected from the model. 

PACS numbers: 61.44.Br, 68.35.Bs, 68.37.Ef, 61.14.Hg 



I. INTRODUCTION 

More then ten years ago, the discovery that centimetre 
size samples of decagonal d-AlCuCo and icosahedral i- 
AlPdMn could be grown opened up the possibility of sur- 
face studies of these quasicrystalsEl. Since then other qua- 
sicrystal samples have been grown to similar dimensions. 
To date, most surface studies haYfj-,hcftjp_n|CiJjMm(^^^ 
the five-fold surface of z-AlPdMnBamOffllMMOBy. 
A consensus has emerged from these studies that this 
surface, after fairly standard ultra-high-vacuum (UHV) 
sputtering and annealing procedures, is itself quasicrys- 
talline. In this work, using a combined experimental and 
theoretical approach, we show that this surface can be 

fnngirlprprl t" hp a tprminatinn of thc kuOWU bulk StrUC- 

ture m^HBEaEl. 
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The dynamical low-energy electroruiiffraction (LEED) 
analysis carried out by Gierer et aZ.E3a indicated that 
the five-fold surface of the t-AlPdMn quasicrystal re- 
tained the bulk quasicrystallinitjtSEi. X-ray photoclec- 
tron diffraction (XPD) studies_ap.|also consistent with a 
quasicrystalline surface natureEJOlld. Large flat terraces 
may be produced, and scanning tunnelling microscopy 
(STM) studies ti3''^6r|jgpseHtfjd similar images pJfr-tii'^ qua- 
sicrystalline surfaceEtmeH-y. Schaub et aimm pro- 
duced detailed STM images of the terraces that reveal 
a dense distribution of dark pentagonal holes of edge- 
length circa 4.8 A oriented parallel to each other, to- 
gether with a more random distribution of bright protru- 
sions. They correlated measurements of structural ele- 
ments both within the terraces and across steps on the 
surface. 

Later, we demonstrated a correspondpice|p£ these mea- 
surements with the geometric model AIESEJO for atomic 
positions of an F-phasecSl. The model M is Jaased on the 
three-dimensional icosahedral tiling T*(^p p]|d<!je«rated 
essentially by Bergman/Mackay polytopesllsOta'cj. The 
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observed terrace structure of the surface was explained in 
terms of the layer structure of the bulk model. The dark 
pentagons observed-on the surface corresponded to the 
Bergman polytopealJ in the bulk layers. The position of 
a given type of terrace was matched to a layer charac- 
terized by a density of certain Bergman polytopes and 
their distribution pattern. We assumed that the surface 
termination respects the integrity of the Bergman poly- 
topes as clusters, at least in the most dense layers, and 
we supposed that such a layer of Bergman polytopes is 
exactly below the termination. However, under these as- 
sumptions it was not possible to explain the observed 
edge-length (circa 4.8 A) of the dark pentagonal holes, 
as this was bigger by the factor r = (VS + l)/2 than 
the pentagonal surfaces of the Bergman polytope (circa 

3 Ayti. 

Later Shen et al., using an autocorrelation analysis 
showed that the surface structure is consistent with a 
bulk structure based on truncated pseudo-Mackay icosa- 
hedraLl and (therefore) Bergman clusters, 
limit of those previous STM studies (RefsE 
), was that the resolution of the images, while sub- 
nanometre, was not atomic. Therefore direct compar- 
ison with bulk models was not straightforward. Addi- 
tionally, the presence of bright protrusions disrupted any 
attempted tiling, and so comparison with tiling models 
was not possible. In a previous paper, we reported an im- 
proved sample preparation technique. This led to a more 
perfect surface devoid of protrusions (section III), and 



in section II B, we analyze in section II D the possibility 



this in turn led to improved resolution in the STM im- 
ages. The better resolution, together with the structural 
perfection, allowed us to demonstrate that the_surface 
structure is consistent with a bulk terminationO, using 
the bulk model of Boudard et a/JlS. 

In this paper we try to find the position of termina- 
tions in the bulk model demanding (i) that the ter- 
minations are ordered in a decorated Fibonacci sequence 
(sections II A, pT^ ) as in RefsoBEl, and (u) that the 
average density of terminations is 0.136 atoms per A^, as 
determined by Gierer et aLEjO (section |[ID| ) . The atom- 
ically resolved images of the surface that allow us t o map 
the local patterns of the STM images (sections [II, [V A ) 
to the local atomic configurations in the terminations in 
M (section IV A ), also prove our Ansatz from section 
[ID which fixes the position of the bulk termination to 
be 4.08 A deeper within LhcJayjcr of Bergman polytopes 



than we expected in With this new position 

of the termination, the edge-lengthj-af the dark pentago- 
nal hole s obse rved by Schaub et a^Jlj is now understood 
(section [V A ) . Moreover we conclude that the termina- 
tion is highly dense in dark pentagonal holes that we can 
now interpret as dissected Bergman (cB) polytopes (sec- 
tion [V A). For the imaged surface terraces the densities 



of single atomic planes in corresponding terminations are 
given and their positions w.r.t. the Bergman layers are 
discussed (section \\ k ) . 

From the general knowledge on possible tilings and 
coverings in five-fold planes in the model developed 



of the existence of the PI tiling on model terminations. 
From the predicted coding windows of t he top five-fold 
q-planes in terminations in M. (section [V A ) we con- 
clude that the Penrose (PI) til ing is possible on almost 



all five-fold terraces. In section [VB we superimpose ex- 
act patches of the PI tiling on STM images of two large 
terraces. 



II. THEORETICAL BACKGROUND 
A. The Geometric Model M 

A geometric model-Jj4, for the atomic positions of i- 
AlPdMn or z-AlCuFe£llll3 has been used ixiinterprct the 
STM measurements data pfjSchaub et alWa on the five- 
fold surfaces of ipAlPdMrO'E3. In the model th£^- 
phasecj (see Ref£3) three-dimensional tiling T*^'^^^t3 is 
decor|a±pd by Bergman (and automatically Mackay) poly- 
topesEllla. For-details on Bergman and Mackay poly- 
topes, see Ref.E^. Thfi— geometric model is based ©a, 
the Katz-Gratias modeler that is explained by Elserll3 
in a three-dimensional "parallel space" , E|| , the space 
in whicij the model projected from Dq latticec3 (see 
RcfsBBl exists. The atoms of i-AlPdMnO or of i- 
AlCuFecj can be placed on three translational classes of 
atomic positions with respect to the Dg-lattice, and are 
denoted qjj^ {= l), b and a, see Table | in section [I C 
and Ref.c2l. These atomic positions in Ey are coded by 
the corresponding "windows" or "acceptance domains" in 
the three-dimensional "perpendicular space", E_l. Note 
that liie six-dimensional Dq lattice, which models an F- 
phaseEj, acts in the six-dimensional space that is a sum 
of E|| and E^. These windows in K± are denoted by Wq, 
Wb and Wa, respectivel}i-|^Hie, windows of the model M 
were constructed in Ref£Zlll3'EJ. The tihng T*'-^^'^ defines 
the quasiperiodic structure. More accurately, the model 
M is supported by tT*^'^^\ the tiling T*(^^) scaled by 
the factor r = (\/5 + l)/2. The quasilattice points of 

All points of the quasilattice which contains the ver- 
tices of the tiling T*^^^) can be embedded in a sequence 
of planes orthogonal to the five-fold symmetry axis of an 
icosahedron ("five- fold direction"). The planes orthog- 
onal to the axis are the "five-fold planes". The planes 
appear in a sequence and have been classified (by partic- 
ular coding regions in the window W7-.2J?) into five types, 
±1, ±2, ±3, ±4, ±5, see Ref.N. The planes of type 1 to 
4 are ordered in the Fibonacci sequence with intervals 
m and I. If the planes of type 5 are included, the se- 
quence of planes forms a "decorated Fibonacci sequence" 
with separations s (short), m (mediu m), I ( long), where 



I — Tm = T s, see Fig. O in section [VA. How is the 



decorated Fibonacci sequence defined? Let us consider 
the Fibonacci sequence of intervals M = I and L = tI. 
If we rename M by / and "decorate" the interval L by 
two points, such that L — m[Js[Jm, the decorated Fi- 
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bonacci sequence appears. For i-AlPdMn that has the 
standard distance parallel to the five-fold direction is ® 

= 4.56 A, and is modeled by tT*^'^^\ s = >' 



2 

T+2 



© = 4.08 A, ; = t2 



2 

T+2 



= 2.52 A, TO = 
= 6.60 A. 

In thejilanes of type 1 a quasiperiodic tiling T*'^'^'')ill 
appearsEEtcEl scaled by a factor t. In the planes of type 
2, 3 and 4, fragments of the sajiie tiling of a plane by 
golden triangles appear (see Ref.llj, Fig.7)|-«fj±ii the same 
inflation properties as in the tiling T*'"^*^ 



In Refdl3^Eil the model is compared to the ideal icosa- 
hedral monograin under the assumption that the terraces 
on the surface of the material are like the planes in the 
bulk, i. e. not reconstructed. This we will first assume 
and then supporijin this paper. The terraces observed 
by Schaub et alBa were related to the sequence of the 
planes of the model Ai described above, see also Ref.E3. 
Whereas Schaub et ai, after annealing at « 800° C ob- 
served onlv-|Eibonacci ordered step heights to and I on 
the surfacOTcl, Shen at ai, after annealing at w 630° C 
also detected the step height s, see Ref.El. 

In this paper we study the fine structure (local atomic 
configurations) within the observed terraces and coippare 
it to the geometric model A^Ej. Whereas in Ref.El we 
succeeded in relating the.|S|equence of the terrace-like five- 
fold surfaces of Schauboa to the layers of the Bergman 
polytopes in the geometric model A^, in this paper, using 
high resolution STM images of a five-fold surface we will 
fix the position of the planes within the layers of Bergman 
polytopes. 

In order to recognize and identify the fine structure 
of the observed surface, we consider certain tilings in 
the five-fold planes and a covering with a set of pro- 
totiles among which are the pentagons and pentago- 
nal starscZI. These tilings will be locally derived from 
the tiling T*^'^^\ The local derivation will be exact 
to a certain stage, and thereafter, random. The tiling 
^*(A4) s(.aled by the factor r defines the quasiperiodic 
structure of the planes on £ke surfaces according to the 
model A4 introduced abovaij. The prototiles in the tiling 
T*i^i) are golden triangles. The edges of the triangles 
in the tiling are parallel to the two-fold symmetry axes 
of an icosahedron ("two-fold directions") and are of two 
lengths, @ and t@. The three-dimensional model Ai is 
supported by the tiling tT*^"^^^ and consequently in the 
five-fold surfaces by rT*^"^*-'. Hence the edges are t(. 



and T @. With the standard value 



4.795 A in the 



case of i-AlPdMn, r@ = 7.758 A and t^(2) = 12.553 A. 

The structure on the surface observed by STM can be 
tiled uniquely only if the tiling, as an abstract structure, 
is derivable from the set of quasilattice points, and if 
the rules of the local derivation are defined on relatively 
small distances with respect to the area of the observed 
surface. 



B. Tilings and Coverings with Pentagonal 
Prototiles contained in the Tiling 7"*^^*' 

As an intermediate step we locally derive the tiling 
T**-^^ with pentagon, acute rhombus and hexagon as pro- 
totiles from the quasilattice T*('^*\ as shown in Fig. [|. 
The tiling has an inflation factor r. It is clear that the 
tiling T**^^) can be reconstructed from its own quasilattice 
points. All edges of the prototiles in T**^^) are of length 
T@. In the geometric model M the prototiles are aug- 
mented by a factor r, so the edge-length is t^(2) = 12.553 
A. All prototiles of T*^^^ are the unions of golden trian- 
gles of the previous tiling T*^^^\ as shown in Fig. 0. If 
we keep that content, the window of the tiling is identi- 
cal to the window of T**^"^"*^ (because none of the vertex 
(quasilattice) points is omitted). The coding window of 
the tiling T*^^\ without content of golden triangles, is 
shown in Fig. ||. Small fractions of the tiling T*^^^ have 
been observed in the five-fold surfaces of decagonal (d)- 
AICuCoEI. 



FIG. 1: The tiling T*'^' of the plane with the acute rhombus, 
pentagon and hexagon as the prototiles. The tiles are marked 
by thick lines and different gray shadows. The tiling 
from which T*'^' is locally derived, is shown in background 
using thin lines. 



FIG. 2: The coding window of the tiling T**^', without the 
content of golden triangles, is inscribed in the decagon by 
thick lines, which is the coding window of the tiling T*^^''-\ 
The codings of the 9 types of vertex configurations in the 
tiling 7"*''**' are marked by the numbers 1-9. 



From the intermediate tiling T**^^^ we can locally de- 
rive a covering of the tiling 7"*(^-i). This covering is by 
two cells in the shape of pentagons, the smaller one, Z3j^, 
of edge-length @ and the bigger, Z3{j, of edge-length t@, 
as shown in Fig. ||(a). Let us denote this covering of the 
tiling T*(^-*) byC^.(^^,. Each acute rhombus from T*^'') 
is transformed into a pair of pentagons of edge-length 
(2) (shown in the left-hand side of Fig. ||(a)), and each 
hexagon is transformed into a pair of overlapping pen- 
tagons of edge-length t@ (right-hand side of Fig. ^(a)). 
The remainder of the tiling T*^^*') should be covered by 
pentagons of edge-length r(2) as in the tiling T*^^\ see 
Fig. ^ The above-defined covering C^,(^^) of the tilin 



tiling 

Mb. 



q-*(A4,) sub-covering of the covering of Kramer 
Kramer also covers the tiling T*^^*' by two pentagons 
of the same size as above. These cells are projected De- 
lone ceUs D°- and of the lattice A4 in E||. In Ey 
they are denoted by ZJj^ and Dy, respectively. Let us de- 
note Kramer's covering by the symbol C^.(A4)- The set 
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of pentagons in C^.(^^) of edge-length r@ is identical 
to the set of Uji's in C^.;^^,. The set of pentagons in 
C^,(Ai) of edge-length @, derived from the acute rhom- 
buses, is a subset of the set of all D'^^'s, in C!^,^a^), and 
therefore the covering C^,(|j,^^ of 7'**^'^*^ is a subcover- 
ing of the covering C^,(^^)E3Eil. Whereas the thickness 
of the covering of C^.,^^, is C'= = 3 - r « 1.382, the 
thickness of the covering of the sub-covering C^.(A4) is 
C = 2t - 2 « 1.236 < 1.382. (For explanation 
of the thickness of the covering see Ref.LS. As a refer- 
ence: a thickness of the covering of a space by a tiling 
always equals 1.) In the subcoverinff C^.^^^, only the 
single and double decking (covering)L3 of the tiles by the 
covering clusters are present. The triple djscking, which 
exists in the covering C^.(A4) is excludecC3 in C^,^^^)- 
The window of the sub-covering C^,;^^) of 7"**^^-*) by 
two pentagons without the content of golden triangles is 
presented in Fig. |[ 



appear obtuse rhombuses. In Fig. ^a) on the right-hand 
side the white spaces around the isolated acute rhom- 
buses are framed by thick lines. Inside these patches, 
there appear pairs of overlapping stars, inscribed in one 
single place in the figure and marked by an arrow. Their 
overlap is exactly the acute rhombus. Up to the choice 
of one star from each pair of overlapping stars, the lo- 
cal derivation of the tiling is exact. The exact tiling of 
the plane by the stars, obtuse rhombuses and pentagons, 
T**^"\ is uniquely determined by its window inscribed in 
the window of T*(^-*\ see Fig. |5|(b). It is the window 
of Niizcki tiling. We randomly choose a star from each 
overlapping pair of stars indicated in the bottom part of 
Fig. H and obtain a partly random tiling T**^"''). The 
only edge-length that appears in the tiling is r@ (in the 
geometric model M, it is = 12.553 A). It is also the 
tiling that could be eventually seen/reconstructed from 
the STM images of the surfaces orthogonal to the five- fold 
direction in i-AlPdMn and i-AlCuFe, and in a decagonal 
phase. 



FIG. 3: (a) The derivation T 



of the figure; (b) T*^^'- 
figure. 



is in the top part 



j^*(pi)i- jg jjj ^jjg ijottQfn part of the 



From the tiling T*^^' let us keep all acute rhombuses, 
and replace each hexagon by two overlapping pentagons 
(as in the sub-covering C^.^a^))- This is an exact lo- 
cal derivation, shown in the left-hand side of Fig. ^(b). 
At this stage we randomly choose one of the pentagons 
from each overlapping pair, and the rest of each hexagon 
unites with the neighboring acute rhombus. In this way, 
either a crown or a pentagonal star appears to replace the 
rhombus, and we obtain a partly random tiling T*^^^)'', 
see the righthand side of the Fig. |^(b) . The ideal class of 
tilings.XPl) with the inflation factor r are described in 
RcfsE^EJ. In Fig. ^(a), the window that exactly defines 
the quasilattice of the tiling T*^^'^'^ is inscribed in the 
window of the tihng T*^^^\ 



FIG. 4: The window of the covering C, 



r'(-*4) 



without the 



content of golden triangles, inscribed by the thick lines in the 
window of the tiling T*'^^^K 



There is another tiling of a plane by pentagonal 
stars^entagons and obtuse rhombuses introduced by Ni- 
izekio. Let us call it the Niizcki star-tiling and denote 
it by T*^"'\ The inflation factor of this tiling is also t. 
In Fig. ^ we derive this tiling from the tiling T*^^\ In 
Fig. ||(a) (top part of Fig. ||) on the left-hand side, from 
the set of all stars, only the locally derivable stars are 
presented. The locally derivable star appears wherever 
there exists an acute rhombus neighboring one or two 
hexagons, each by an edge. Between these stars, there 



FIG. 5; (a) The window of T*'''^' is inscribed in the window 
of 7"*''^*^ by thick lines, (b) The exact tiling of the plane by 
the stars, obtuse rhombuses and pentagons, T*^"\ is uniquely 
determined by its window inscribed by thick lines in the win- 
dow of T*'^*'. It is the window that codes the Niizeki tiling 



Both exact tilings, T*(pi) and T*("), can be locally 
derived from their respective quasilattice points. In the 
reconstruction of tilings T*^^^' and T*^") from the re- 
spective quasilattices there appear: (i) pairs of pentag- 
onal sets of points centered in each other and mutually 
rotated by 27r/10. The set of points of the smaller size 
(the smallest pentagonal set in the tiling) is on the neigh- 
boring distances (2), the bigger, on neighboring distances 
T^@. Each pair leads to the pentagonal star; (ii) the 
isolated pentagonal sets with neighboring distances r@ 
are to be connected in pentagons. In order to reconstruct 
the tiling T*'^^^ , it is enough to draw the pentagons from 
the isolated five-tuples of five-fold symmetrically ordered 
points. In order to reconstruct the tiling T"**^"^ one draws 
the stars from the pairs of pentagonal sets defined above. 
One can show that in an abstract sense the tilings J^p^") 
and T"**^"' can be mapped one-to-one to each otherEa. If 
we consider an experimental atomically resolved five-fold 
surface, and tile the observed surface, we first have to 
identify the surface by a plane in the model which we 
will call an Al-plane. Then we determine the coding 
window of the plane in Ej_, the Al-plane window, and 
we place the biggest possible window of an exact tilings 
7-*(pi) or T**^"^ in the A^-plane window, see Figs. and 
p^ . Following these arguments, we will determine the 
edge-length of a possible tiling of an observed surface by 
the prototiles of the tiling PI in section IV. 
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FIG. 6: Local derivation: T**^) — > T**"' 
top part of the Figure is referred to as ( 
part as (b). 



In the text, the 
and the bottom 



Atomic Positions in Five-Fold Planes of the 
Geometric Model M 



In section II B we have derived the tilings T*(^\ T*^p^^ 
and T*(") either from the ideal tiling T*^"^*^ or from 
their own corresponding quasilattices. We have been 
considering exclusively these points q € Dq that belong 
to the underlying tiling of the model, rT*'"^*"^. Conse- 
quently the edge-lengths in both locally derived tilings, 
r*(pi)'- and T*("-) were of length t^® = 12.553A. If 
we also take into account the decoration of the tiling by 
Bergman/Mackay polytopes, the window of the quasilat- 
tice points of tvpe oG Dq, Wq^^ becomes the polytope 
derived in ReffOli'S. 

In order to study the five- fold planes of the model M , 
we present two important general facts that we implicitly 
use in all our considerations. 

(i) The reciprocal lattice of the root-lattice Dq we de- 
note by Dq^'^'^. The lattice Dq^'^'^ is also known as the 
weight-lattice Dq" . If one icosahedrally projects Dq^'^'^ 
to the parallel ^m;e, E||/ (E_l), an icosahedral Z(t)- 
module appearsHEj. The module points in a plane of a 
three dimensional (icosahedral) Z(T)-module in E||, un- 
der the * - mapE3, i. e. r — > ^'I'e mapped in 

into a plane too. The section of this plane in Ej^ 
through the three-dimensional window (acceptance re- 
gion of the three-dimensional quasilattice) defines a two- 
dimensional window of the quasilattice in a correspond- 
ing plane in E|| . The analogous statement holds true for 
the lines. These are the general properties of a Z(A) mod- 
ule with quadratic irrationality A. In our considerations 
A = r = (-^5 + 1)/2. The above statement is valid for the 
modules with symmetries such as icosahedral, five-fold, 
ten-fold, eight-fold and twelve-fold. 

(ii) Let us consider the four translational classes with 
respect to the root-lattice Dq of six-dimensional points 
i(ni, ...jTig) G Dq, where Ui are integers. The condition 
for an atomic position x = ■^{ni, ...jTig) to be in a five- 
fold plane in E|| or E_l is a class-function presented in 
Table |. Hence, the atomic positions in a five-fold plane 
of a De^-icosahedrally projected Z(r)-module belong to 
the single class, q_Dg(= q),b, a or c. 

Using the facts (i) and (ii) , in the geometric model Ai 
we code each five-fold plane containing a class of atomic 
positions in E|| by the five-fold dissection in Ej_ of the 
single window Wq, Wb or Wa, corresponding to that class. 



TABLE L The condition for atomic positions x = |(ni , ng) 
to be in a five-fold plane in Ey or Ex is a class-function. The 
symbols e and o stand for even and odd integers, respectively. 
The symbol ny / is a unit normal to the five-fold plane 
in E||/ Ex space, respectively, G Ey and x± £ Ex, where 
X is the point in six-dimensional space, Ey-I- Ex. The scalar 
product is given in the units [k], k — l/[\/2(r + 2)]. 
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D. Densities of Five-Fold Planes and Terminations 
of Five-Fold Surfaces in the Geometric Model A4 

In the work of Gierer et an average density of 

"terminations" of five-fold surfaces has been determined 
to be pq+b = Pq + Pb = 0.136 atoms per A^. By the 
density of the termination the authors mean the sum of 
the densities of two atomic planes on a surface separated 
by a vertical distance of 0.48A, and consequently each 
"termination" corresponds to a pair of planes separated 
by this distance. Let us suppose that the surface (top) 
planes are of type q, then the planes 0.48 A below, in 
the geometric model M, are of type b. Let us calculate 
Pq+b(z||) in the model, where zy is along a five-fold axis 
orthogonal to terraces on the surface, and let us plot this 
value along the corresponding z±, Pq+b{z±). The result 
is shown in Fig. |^. 



FIG. 7: Density Pq+b of the pairs of five-fold planes in bulk 
model A4: a q-plane and a b-plane, 0.48A below the q-plane. 
pq+b as a function of zj_ in units of r^ ©. The image of z-axes 
in Ex, z±, is chosen such that zy points into opposite direc- 
tion of the bulk. pq{z±) is the density of a q-plane, pb{z±) 



is the density of a b-plane shifted by c 



Vir + 2)]®, 



Pq{z±) + pb(z±) = pq+b(z±)- In the figure the oWand the new 
coding regions of the (decorated) Fibonacci sequence of planes 
that represent the surface terraces in M are marked. In the 
new region, th« representative plane of the biggest clear ter- 
race of Schaubu is marked by S8 on a new position. The condi- 
tion for appearance of the Ring-plane in a sequence ml(R)ml 
is determined and a representative of a Ring-plane (R) to- 
gether with a representative of the following Clear-plane (C) 
are marked. Finally the region of existence for PI tilings on a 
q-plane is denoted and particular minimal edges are attached 
to their coding regions, see section IV, 



The function Pq+b{zi_) has clear (almost fiat) plateau. 
The appearance of the plateau is due to the polytopal 
shape of the codup^hpAandows Wq and Wb in the geo- 
metric model AIIIJ'ESEJ. In particular the window Wq, 
that defines the surface (top) plane in a termination, dif- 
fers strongly from the spherical shape. The plateau of 
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the function pq^b{z±) simultaneously contains the max- 
ima of the function and has a value that approximately 
equals the average density of terminations determined 
by Gierer et a/.O, which is 0.136 atoms per A^. It is 
easy to conclude that all terminations on terraces must 
have equal densities. Consequently an interval on z±_ un- 
der the plateau (the "carrier" of the plateau ) codes the 
terminations and indicates the new coding of surface q- 
planes |tDjiie_shifted from the pZcLjialiie that we expected 
in RefsooEil. In those papersoo'Eil we supposed that a 
single surface q-plane has to have the highest density. In 
accordance with this Ansatz, at least some of the dense 
layers of Bergman poiytopes were below the. surfaces. 
But the dark pentagons observed by Schauta, that we 
put in correspondence to the Bergman poiytopes in the 
layer below the surface, were_higger by a factor t than the 
faces of Bergman polytopealj. Let us shift the surfaces 
of terraces by 4.08A in the direction of the bulk (-4.08A 
along Z|| ) in "parallel" (observable) space, such that the 
q-plane on a terrace dissects the Bergman poiytopes of 
the layer and the section of each Bergman polytope is a 
pentagon of edge-length 4.8A, approximately the size of 
the dark pentagons observed by Schaubcl. This shift in 
E|| corresponds to the shift by [2t/ (r -I- 2)]® along z± in 
orthogonal space. Indeed, the coding interval of q-planes 
that forces the planes on the surface to appear in a Fi- 
bonacci sequence (or in a decorated Fibonacci sequence) 
is placed under the plateau of the function /9g+h(z^) by 
this shift, see Fig. |^. We suppose that the terrace-like 
five-fold terminations do appear in a Fibonacci (or deco- 
rated Fibonacci) sequence such that the top q-planes in 
terminations need not be the most dense among the q- 
planes, but the above defined "terminations" , the pairs of 
planes on a surface, have the highest densities among all 
such pairs of q- and b-planes in the geometric model M . 
We check our hypothesis (Ansatz) on two large terraces 
in section 



angle followed by annealing for periods of about 1 hour 
at 970 K. The results are shown in Fig. ||(left hand panels 

(a) and (c)). For the second preparation. Preparation II, 
a further polish using 0.25 /im diamond paste was used. 
The surface was prepared in- vacuum by several cycles of 
sputtering with 0.5 keV Ar ions, with a sputtering angle 
of 20°-30° relative to the surface parallel, followed by an- 
nealing to 970 K for two hours (in total twelve hours of 
annealing); Fig. ^(right hand panels (b) and (d)) show 
the results. 

When large scale scans are compared (Fig. ^(a) and 

(b) ), it is evident that larger terraces are obtained using 
Preparation II. For Preparation I, the largest terraces 
are of the order of 1200 A in magnitude. For Prepara- 
tion II terraces of width 4000 A and length of micron size 
were obtained. Further differences between the results 
of the preparation techniques are observed when scans 
of smaller area are compared. Fig. ||(c) and (d) show 
100 A X 100 A areas of each surface. Clearly the surface 
in Fig. ||(c) is not as well resolved as that in Fig. ^(d); 
the bright spots in Fig. ||(c) correspond to protrusions of 
height up to 2.0 A, while dark spots are associated with 
holes of depth estimated to be at least 1.5 A. This STM 
im*ise,|ifi comparable to those in the work of Schaub et 
al.auuQ. This can be contrasted with the surface shown 
in Fig. |^(d) where there are no large protrusions and the 
surface corrugation within the terraces is < 1 A. Because 
the STM tip can scan the surface more closely, features 
on the surface are better resolved. The features in this 
image have dimensions typical of atomic sizes (2-3 A). 
Larger features (4-6 A) are also evident and probably 
represent groups of a few atoms. The LEED patterns 
from each of these surfaces are qualitatively identical, 
but the range of electron beam energies over which the 
LEED patterns are obtained is much larger (10-300 eV) 
using Preparation II than for Preparation I (40-180 eV). 
The LEED patterns have very sharp diffraction spots, a 
low background, and show five-fold symmetry. 



III. FIVE-FOLD SURFACES IMAGED BY STM; 
SURFACE PREPARATION AND STM 
RESOLUTION 

In this section we describe the surface preparation we 
have developed to obtain large flat terraces and low sur- 
face corrugation in STM experiments. We contrast STM 
results using our optimum preparation wilj^presults pre- 
viously published by us and other groupsaLru. 

Fig. H shows data from the surface of i-Al-Pd-Mn af- 
ter the two different preparation procedures. In each 
case the quasicrystal samples were gr own , at Ames Lab- 
oratory using the Bridgman methodEllEa. After being 
cut perpendicular to their five- fold SMnmetry axes in air, 
the sample surfaces were preparedlj by polishing. For 
the first preparation. Preparation I, the sample was pol- 
ished using 6, and 1 fim diamond paste for one hour. 
In- vacuum preparation consisted of a few cycles of argon 
ion sputtering at 1 keV energy and a normal incidence 



FIG. 8; (a) 1500 A x 1500 A STM image showing atom- 
ically flat terraces from a surface prepared using Prepara- 
tion I. (b) 17500 A X 17500 A STM image showing atom- 
ically flat terraces from a surface prepared using Prepara- 
tion II. (c) 100 A X 100 A STM image of a flat terrace that 
we call the "Clear", C-terrace from a surface prepared using 
Preparation I (bias voltage 2.29 V, tip current 0.59 nA). (d) 
100 A X 100 A high resolution STM image of the same C- 
terrace obtained on the five-fold surface using Preparation II 
(V=l V, 1= 0.3 uA). 



The resolution can be put on a semi-quantitative basis 
by calculating the two-dimensional lateral autocorrela- 
tion functions of the images of Fig. ^(c) and (d). These 
are shown in Fig. ^(a) and (b) respectively. While the 
symmetry of both autocorrelation patterns is similar, the 
pattern of Fig. ^(b) is considerably clearer and the cor- 
relation maxima extend to longer distances indicating a 
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higher degree of quasiperiodic order. 

For a more quantitative comparison a radial distribu- 
tion function (RDF) has been calculated in both cases. 
The procedure consists of dividing the 360° around the 
centre of the autocorrelation function in increments. 
Along each line corresponding to each increment, the dis- 
tances from the centre to the maxima are measured. All 
the measurements are then averaged and plotted as his- 
tograms (Fig. ^(c) and (d)). It can be seen that there 
is considerably more structure in the RDF in Fig. ^(d) 
than in that of Fig. ^(c). 



FIG. 9: (a) 100 A x 100 A lateral autocorrelation function of 
the STM image of Fig. |(c). (b) 100 A x 100 A lateral auto- 
correlation function of the STM image of Fig. ^(d) . (c) Radial 
distribution function calculated from the autocorrelation pat- 
tern of (a), (d) Radial distribution function calculated from 
the autocorrelation pattern of (b). 



In summary, surfaces prepared using Preparation II 
have a much lower surface corrugation and lead to much 
better resolved STM data than those previously obtained 
using Preparation I. The main differences in these proce- 
dures are the sputtering energy and incidence angle (sug- 
gesting that minimizing surface damage while removing 
contaminants is of importance) and the long anneal times 
at high temperatures which probably serve to restore the 
surface composition to that of the bulk quasicrystal. We 
interpret the protrusions as due to material on the surface 
which has not yet diffused to the step edges. A similar 
observation was recently made to explain the origin of 
such protrusions on d-Al-Ni-Co surfacestj. 



IV. REPRESENTATIONS OF SURFACES ON 
TERRACES IN THE GEOMETRIC MODEL M; 
TILING ANALYSIS OF STM IMAGES 

A. A Five-Fold Terraces Mapped to the 
Terminations in M 



In section II D we suggested new positions of five-fold 
terminations in the geometric model A4. In this section 
we search for the terminations in the geometric model 
^A on these new positions that fit to atomically resolved 
pictures of five- fold surfaces on particular terraces imaged 
by STM. In a sequence of five-fold terraces we observe a 
large terrace that contains a rare local configuration that 
we call the "Ring" (R). This configuration helps us to 
orient in the bulk model A4, i. e. to fix the position of 
the R-terrace w.r.t. the five-fold z-axes. Near the R- 
terrace we observe the clearest terrace that we denote by 
"C". A fragment of the C-terrace is shown in Fig. ||(d). 

On the C-terrace local configurations of the five-fold 
depressions in the shape of dark Stars (dS) are observed. 
The strongly shining pentagonal local configurations in 



the form of the white Flower (wF) and the white Star 
pointing upwards (wSu), both parallel to the dS and in 
the same direction make a white picture on a dark back- 
ground, see Fig. [lo|(a). 



FIG. 10; (a) 100 A x 100 A high resolution STM image of the 
C-terrace on a five-fold surface. On the C-terrace frequently 
repeated local configurations such as a dark Star (dS), a white 
Flower (wF) and a white Star pointing upwards (wSu) par- 
allel to the dS are marked. The Bergman polytope below 
the terrace (Bb), above the terrace (Ba), and the Bergman 
polytope dissected by the terrace (cB) are also marked. For 
the scale the wSu is framed by a pentagon of edge-length 
tD, D ~ 4.8A. (b) The C-terrace from (a) corresponds to 
the C-termination in M. Black points are atomic positions 
in the q-1024 plane in M (No 175 on Fig. |l^) which is on 
the surface, grey points are in the b-1025 plane, 0.48A below 
the q-1024 plane. The local atomic configurations that may 
present the dS, the wF and the wSu are marked. The main 
constituents of these configurations are the top surface of the 
Bergman polytopes that are in the layer below the surface 
(Bb), the bottom surface of the Bergman polytopes that are 
in the layer above the surface (Ba) and the pentagonal section 
of the Bergman polytopes from the layer that is dissected by 
the surface (cB). Scale: D = rd — 4.8A. From the center of 
cB the next atomic position in the bulk is 2.04 A below the 
surface. 



In contrast to the C-terrace the R-terrace is not con- 
tinuously (globally) clear, i. e. the STM images of the R- 
terrace taken on different places lead to different RDFs. 
Nevertheless, we observe some local configurations on the 
R-terrace that are clear, see Fig. |ll|(a). We find the white 
Flower (wF) and the dark Star (dS) identical to those on 
the C-terrace, see Fig. |l^(a), but we also see a charac- 
teristic "Ring" -configuration (R) that is present on none 
of the other observed terraces. The terrace is therefore 
denoted the R-terrace. In the R-terrace there is also a 
configuration which we call the white Star pointing down- 
wards (wSd), it is rotated 180° w.r.t. the wSu that we 
observe on the C-terrace. 

As we stated the areas of both R- and C-terrace are 
large and they appear in a local upward sequence of 
steps ml(R)ml(C), where m « 4. 08 A and I « rm. On 
the q-planes of the geometric model M we find a rare 
atomic configuration that may represent a local Ring- 
configuration on the STM-image of the R-terrace, com- 
pare Figs |l^(a) and (b). We determine the coding of 
the ring-configuration (R) in and demand that the 
q-plane containing the Ring-configuration is to be found 
in an upward sequence of the q-planes ml(R)ml( C) o n 
the new positions (shifted by 4.08A, see section \1 D| ), 
and both R- and C-plane are to be among the planes 
from the decorated Fibonacci sequence. From these con- 
ditions we find the coding area in along z_\_ of the 
R-plane to be in the interval e (0.198, 0.337)[t2®] 
marked on Fig. |^. In a patch of the geometric model 
that spreads along zii-axes in an interval of 1195 A 
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FIG. 11: (a) 75 A X 75 A STM image of the R-terrace. The 
local configurations Ring (R), dark Star (dS) and white Star 
pointing downwards (wSd) are framed by three pentagons of 
edge-lengths t^D, t^D and rD respectively, where D ~ 4.8A. 
On a bigger STM image of the R-termination a full white 
Flower (wF) can be seen also, (b) the R-terrace from (a) 
corresponds to the R-termination in M. Black points are 
atomic positions in q-1037 plane in M (No 178 on Fi| 
which is on the surface, grey points are in b-1038 plane, 
below the q-1037 plane. The local configurations of atomic 
positions that may represent the dark Star (dS), the white 
Flower (wF) and the white star (wSd) anti-parallel to the 
dark Star (dS) are marked. Scale: D = rd = 4.8A. In the 
centre of the dark Star the nearest atomic position is 2.04 
A below the surface. In the q-1037 plane there are empty 
"streets", A — 4.56A broad. 



we find only 15 representatives of the R-plane that fulfiU 
all conditions mentioned above. We choose the q-1037 
plane (No 178 on Fig. |l^), the plane that is coded in 
Ex by z± — 0.323t^®, see Fig. ^ The corresponding 
C-plane is then q-1024 (No 175 on Fig. ^ coded by 
z± = 0.192r^@, see Fig. 0. In Fig. |l2| the coding win- 
dows of the R- and the C-plane are shown. This pair 
of R- and C-planes (one of 15 pairs in the model-patch) 
are taken not far from \he estimated model-plane S8 for 
Schaub's terrace No. 8lI ( 



on a new position q-1128 (No 
193 on Fig. |13|) coded by z± = 0.211r2@, see Fig. g 



FIG. 12: In Ex the windows of the top (q-)planes in R- 
and C-terminations, Wr and Wc, respectively. Over them is 
plotted (i) the window of the tiling PI of edge-length 4.8A (in 
E||) denoted by W/^pi). It is the maximal window of PI, such 
that VF(pi) C Wc and (ii) the window of the tiling r(Pl), 
of edge-length 7.8A (in Ey), denoted by Wt-(pi). It is the 
maximal window of PI, such that Wtpi C Wr{C Wc)- The 



scale for the figure is set by the decagon W, 



which is 



the window of the tiling T*'-'^*'-' with edges d = t~^D and 
D = 4.8A (in Ell) . For the biggest possible window of PI in 
W^,(A4) see Fig n(a) in section II B. 



FIG. 13: A decorated Fibonacci sequence (s = 2.52A, m = 
4.08A, I = 6.6OA) of the q-planes along the (five-fold axes) 
of t|fj)e ±1, ±2, ±3, ±4, ±5 in M on the old positions, see 
Ref.lij. Relative to these positions the stacked layers of the 
Bergman polytopes are drawn with their relative densities. 
The representative planes of the large R-, C- and S8-terraces 
are marked on the new positions. The — 4.08A shift from the 
old to the new positions is indicated by arrows. 



As we have stated, in contrast to the R-terrace, the C- 
terrace is uniformly clear and it has a unique radial distri- 
bution function (RDF). Fig. 11^ (top) corresponds to the 



FIG. 14: Radial distributions calculated (top) from the auto- 
correlation pattern of the high resolution STM image shown 
on Fig. pi d) ; (bottom) from the autocorrelation pattern of the 
q-plane (q-1024) of the C-termination in the geometric model 
A4 presented in Fig. p^(b) 



RDF calculated from the high resolution STM image of 
Fig. |(d) (identical to the RDF shown in Fig. ||(d)). Max- 
ima are found at 7.3, 12.1, 19.4, 24.2, 31.1, and 38.0 A 
(± 0.3 A). The radial distribution function calculated 
from the C-plane q-1024 (No 175) of the geometric model 
Ai (shown in Fig. [T4|(bottom)), is very similar, the main 
differences being the presence of a double peak at 15 A, 
and some extra structure at higher distances. The cor- 
respondence with the largest intensity peaks is however 
very good. 

To the C- and the R-planes of type q, there corre- 
spond C- and R-terminations, which are pairs of q- and 
b-planes at the surface separated by a vertical distance of 
0.48A. All local patterns observed on the C-terrace and 
the R-terrace can be mapped to the model-terminations, 
the C(lear)-termination and the R(ing)-termination, re- 
spectively, see Figs |l^(b) and |ll|(b). These patterns 
are mapped to the local atomic configurations on model- 
terminations that contain groups of atoms in the shape 
of pentagons related to Bergman polytopes that (i) ei- 
ther dissect by the termination (cB) (some of them are 
in the central parts of dS); or (ii) are below the termina- 
tion (Bb) (some of them are in dS, wF, R, wSd); or (iii) 
to the Bergman polytopes that are above the termination 
(Ba) (in wSu). For this see Figs |l^(b), |l^(b) and |l|. 



B. The Tiling PI on Five-Fold Surfaces 

In order to extract information from thC[-S3jM images, 
we have employed a tiling approach in Refa3'c2l. In Ref .c3 
this consisted of connecting points of high contrast on 
the STM image to create pentagons. The fiUing-in of the 
image using pentagons led to a Penrose (Pl)-like tiling 
of the experimental plane (with an edge-length of 7.8A). 
Here we will reconstruct exact patches of the PI tiling on 
the STM images of both R- and C-terraces (see Figs n3 
and |l6|) and on corresponding model-planes (not shown]^^ 

The coding regions of PI tilings with minimal edge- 
lengths on q-planes are marked on Fig. ^. The tiling 
PI with edge-length 7.8A is coded in the interval zj_ e 
{~t^^,t^^)[t'^©] and that with edge-length 4.8A in the 
interval z± e (-r"^, r~3)[T^®]. 

From the coding of the q-planes of the C-, S8- and 
R-terminations {z^ = 0.192r2®, zfs = 0.211r2® and 
= 0.323t2®) we conclude that the q-1024 plane (No. 
175 in Fig. |l^) of the C-termination and the q-1128 plane 
(No. 193 in Fig. |l|) of the S8-termination in M allow 
a PI tiling of minimal edge-length 4.8A, and the q-1037 
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plane (No. 178 in Fig. |T^) of the R-termination allows a 
PI tiling of minimal edge-length 7.8A. (See also in Fig. |l3 
the coding windows of PI tilings with edge-lengths 4.8A 
and 7.8A plotted over the coding windows of the q-1024 
plane of the C-termination and the q-1037 plane of the 
R-termination) . 

An exact patch of the tiling PI can be exactly placed 
on the q-plane of a model termination as follows: (i) Plot 
the window of the PI tiling, Wpi, of the maximal possible 
size such that Wpi C Wq-pi, where Wq-pi is the coding 
window of the surface q-plane in the model M . For the 
biggest possible window of PI in Wq-.iA^-^ see Fig ||(a) 
in section II B. (ii) Mark all atomic positions coded by 
the points in the window Wpi in E|| . This set of points 
uniquely determines the PI tiling on the model q-plane. 
The procedure is evident. In contrast, for an STM image 
of a terrace we have to proceed locally. If the plane is 
very clear and the window Wpi can be tightly placed in 
the corresponding window Wq-pi we can reconstruct an 
exact patch of the PI tiling by trial and error. A probable 
exact patch of the tiling PI with minimal edge- length of 
7.8A is reconstructed on an STM image of the R-terrace, 
in Fig. M. 



FIG. 15: 75 A X 75 A segment of an STM image of the 
R-terrace with a superimposed exact patch of the PI tiling of 
edge-length 7.8A. 



FIG. 16: (a) An exact patch of the PI tiling superimposed on 
the enhanced high resolution STM image (100 A x 100 A) of 
the C-terrace. (b) The patch of PI tiling obtained from (a) 
shown superimposed on the unenhanced high resolution STM 
image of the C-termination from Fig.^d). 



to strongly enhance features in the image correspond- 
ing to the selected fc- values. The procedure is essentially 
identical to that used by Beeli and co-workers in the en- 
hancement of High Resolutkffl Transmission Electron Mi- 
croscopy (HRTEM) imagCEcil. 

In the enhanced image the white spots that we inter- 
pret as the images of atomic positions are almost as sharp 
as in the model plane q-1024 from the C-termination, see 
Fig |o|(b). We find a patch of exact PI tiling of edge- 
length 7.8 ± 0.2 A that can be easily superimposed on 
the enhanced image, see Fig. |l|(a). Fig. |l^(b) shows 
this tiling superimposed on the unenhanced STM image. 



C. Densities of Five-Fold Planes and Terminations 

In Table |l| we compare the densities of the R-, C- 
and S8-terminations, and also the densities of single (q- 
and b-) planes contained in each termination on the old 
and the new positions. It is evident that the densities of 



TABLE II: Densities on old and new positions of the R-, S8- 
anri-iC-terminations (shifted by —4.08 A). Following Gierer et 
a/.cj a termination contains two planes on the s urface , and in 
M its density is P{q+b) ~ P(q) + P(6), see section II D, 



termination 


R 


S8 


C 


average 


No(PI°''') 


177 


192 


174 






-0.019 


-0.131 


-0.150 




^old r A — 2i 


0.087 
0.026 


0.084 
0.008 


0.082 
0.007 






0.113 


0.092 


0.089 


0.098 


No(pi:;™) 


178 


193 


175 






0.323 


0.211 


0.192 




P(,Ti>)[A"'] 


0.059 
0.076 
0.135 


0.074 
0.063 
0.136 


0.076 
0.060 
0.136 


0.136 



The q-1024 plane related to the surface of the C- 
termination is very dense, and although we could theoret- 
ically place the PI tiling of minimal edge-length 4.8A (see 
Fig. |2|), we have managed to reconstruct only an exact 
patch of PI tiling of edge-length 7. 8 A on the STM image 
of the C-terrace, see Fig. |^ For this purpose we apply 
an image enhancement technique to the data of Fig. H(d) 
in order to even out experimental contrast variations (in- 
herent in the use of the STM technique which measures 
electron charge density at the surface rather than nuclear 
coordinates) and to reduce experimental noise. The pro- 
cedure is based on Fourier filtering and consists of taking 
a fast Fourier transform of the image, and then enhanc- 
ing obvious Bragg reflections with unique fc-values and 
removing experimentally-induced diffuse features due to 
noise. This modified frequency space representation is 
then Fourier transformed to obtain the filtered image 
shown in Fig. |l^(a). The result of this procedure is 



terminations on the new positions give a better fit to the 
LEED result of Gierer et aZJlJ, an average density of 0.136 
atoms per A2. We see that the C- and S8-terminations 
contain top q-planes that are much more dense compared 
to the top q-plane of the R-termination. The STM im- 
ages of the C- and S8-terminations show that they are 
continuously clear. Another fact is that the top q-plane 
of both (C- and S8-) terminations are on similar relative 
positions in the bulk w.r.t. the layers of Bergman poly- 
topes (compare Fig. |l^): a dense layer (Bl) is below the 
plane, a middle dense layer (B') is dissected by the plane 
and a layer of low density (B") is above the plane. 



V. CONCLUSIONS 

We have presented two atomically resolved, high reso- 
lution STM images of large and flat terraces on the five- 
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fold Al7oPd2iMn9 surface. We have mapped these sur- 
faces to the five-fold terminations in the geometric model 
M such that they form a decorated Fibonacci sequence, 
and their average atomic density is iit^greement with the 
LEED measurements of Gierer et alS^. Due to the poly- 
topal windows of the geometric model A4 all terminations 
turn out to have equal and simultaneously maximal den- 
sities. These new terminatiotus in Ai are placed 4.08 A 
lower than in the work of Ref.tj. On the present STM im- 
ages the dark pentagons appear as the dark Stars. At the 
new positions of the model-termination planes the pcjt-i 
terns of dark pentagonal holes are the same as in Ref.t2l 
but now each dark hole is of an appropriate size. At 
the new positions the surface terminations dissect the 
most dense Bergman layers in the model Ai. The local 
patterns on STM images are present in the model termi- 
nations and are related to the Bergman layers above (if 
one exists), below and dissected by the termination. Dis- 
sected Bergman polytopes correspond to the dark Stars. 
The edge-lengths of superimposed exact patches of the 
Penrose PI tiling on two STM images (corresponding to 
pentagons of height equal to 12 ± 0.36 A) arc shown to 
be in agreement with the bulk model A4 based on the 
tihng tT^^^M. 



VI. ADDENDUM 

Since this manuscript was submitted for publication 
another paper containing STM results has been pub- 



lished, see Refii. We note that those authors also con- 
clude that the i-Al-Pd-Mn surface is a termination of the 
bulk structure. 
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